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ABSTRACT 

We utilize the newly- acquired, ultra-deep WFC3/IR observations over the HUDF to search for star- 
forming galaxies at 2; ~ 8 — 8.5, only 600 million years from recombination, using a Yios-dropout selec- 
tion. The new 4.7 arcmin^ WFC3/IR observations reach to ~28.8 AB mag (5cr) in the 1^105^125^160 
bands. These remarkable data reach ~1.5 AB mag deeper than the previous data over the HUDF, 
and now are an excellent match to the HUDF optical ACS data. For our search criteria, we use a two- 
color Lyman-Break selection technique to identify 2: ~ 8 — 8.5 yio5-dropouts. We find 5 likely z~8-8.5 
candidates. The sources have i^ieo-band magnitudes of ~28.3 AB mag and very blue [/F-continuum 
slopes, with a median estimated of < —2.5 (where fx oc A^). This suggests that z ^ 8 galaxies are 
not only essentially dust free but also may have very young ages or low metallicities. The observed 
number of Fio5-dropout candidates is smaller than the 20±6 sources expected assuming no evolution 
from z ~ 6, but is consistent with the 5 expected extrapolating the Bouwens et al. 2008 LF results to 
z ^ 8. These results provide evidence that the evolution in the LF seen from z ^ 7 to z ^ 3 continues 
to 2; ~ 8. The remarkable improvement in the sensitivity of WFC3/IR has enabled HST to cross a 
threshold, revealing star- forming galaxies at z~8-9. 
Subject headings: galaxies: evolution — galaxies: high-redshift 



1. INTRODUCTION 

An important uncharted frontier is understanding how 
galaxies build up and evolve from the earliest times. 
While great progress has been made in characterizing the 
galaxy population at z < 6, extending these studies to 
z > 7 has proven extraordinarily challenging. Only ~25 
high-quality z ^ 7 candidates are known (e.g., Bouwens 
et al. 2008; Oesch et al. 2009; Ouchi et al. 2009; Castel- 
lano et al. 2009; Gonzalez et al. 2009; R.J. Bouwens et al. 
2009, in prep). Fundamentally, the challenge has been to 
obtain extremely deep observations in the near-IR, where 
the redshifted UV light of faint z ^ 7 galaxies is found. 

Now, with the installation of the WFC3/IR camera on 
HST, we have a far superior surveying instrument, with 
6x the area of NICMOS, >2x the resolution, and 2-4x 
the sensitivity. These capabilities allow us to search for 
2: > 7 galaxies ~40x more efficiently. 

Here we report on our use of the early WFC3/IR ob- 
servations over the HUDF to search for galaxies at 2; > 8. 
This is the same epoch in which a GRB was recently dis- 
covered at z ~ 8.2 (e.g., Salvaterra et al. 2009; Tanvir 
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et al. 2009). Throughout this work, we quote results in 
terms of the luminosity L^^^ Steidel et al. (1999) de- 
rived at z ~ 3, i.e.. Mi 700 ab = —21.07. We refer to the 
F606W, F775W, F850LP,' F105W, F125W, and F160W 
bands on HST as FeoG, ^775, ^850, ^105, -^125, and i^ieo, 
respectively, for simplicity. Where necessary, we assume 
Qo = 0.3, Qa = 0.7, Ho = 70km/s/Mpc. Ah magnitudes 
are in the AB system (Oke & Gunn 1983). 

2. HUDF WFC3/IR OBSERVATIONS 

The present high redshift galaxy searches utilize the 
first epoch of ultra-deep near-IR WFC3/IR observations 
acquired over the HUDF (Beckwith et al. 2006) for the 
192-orbit HUDF09 program (G011563). This program 
will create three ultra-deep WFC3/IR fields, one po- 
sitioned over the HUDF and the other two over the 
HUDF05 fields (Oesch et al. 2007), each imaged in three 
near-IR bands yio5, ^125, and i^ieo- Combining these 
ultra-deep near-IR data with the similarly deep optical 
HUDF data permits us to select 2; ~ 7 2^850, ^ ~ 8 F105, 
and even z ^ 10 J dropout galaxies to very low lumi- 
nosities (i.e., -18 AB mag, -0.06 L^^s). 

The WFC3/IR field over the HUDF is centered on 
3^32^38.5^ and -27^4^0.0^'. In the first year of ob- 
servations, we obtained 16 orbits of Fio5-band data (2 
orbits were severely impacted by persistence and are not 
included), 16 orbits of Ji25-band data, and 28 orbits 
of i^i6o-band data. The 60-orbit observations were ob- 
tained from August 26, 2009 to September 6, 2009. 

Standard techniques were used to reduce the HUDF09 
WFC3/IR imaging data. Individual images - after mask- 
ing out sources - were median stacked to create super 
median images (one per filter) and these median images 
were then subtracted from the individual frames. The 
images were then aligned and drizzled onto the same grid 
as the vl.O HUDF ACS data (Beckwith et al. 2006) re- 



2 




M25 



160'' AB 



Fig. 1. — Yios — J125 vs. J125 — i^ieo two-color selection used 
for identifying z ~8 lios -dropout galaxies in our HUDF WFC3/IR 
data. The colors of our 5 F-dropout candidates are indicated with 
the large red squares (with Icr error bars and 2a limits: to be 
conservative) and open red squares (with Icr limits). The blue lines 
show the colors expected for young star-forming galaxies (with UV 
continuum slopes ^5 = — 3 or ^5 = —2) as a function of redshift. 
The solid, dashed, and dotted red lines show the colors we would 
expect for various low-redshift galaxy SEDs (Coleman et al. 1980). 
The hatched green region indicates the region in color space we 
would expect L,T dwarfs to lie (Knapp et al. 2004). The colors 
of individual sources in our photometric sample are given with 
the small black points (for those detected in the HUDF optical 
data and hence not included in our selection) and the small red 
points (for those undetected in the optical). The inset to the figure 
presents the expected redshift distribution for the sample using the 
simulations described in §4.1 and the fiducial LFs from §4.2. All 
of our z ~ 8 Yios-dropout candidates have bluer t/y-continuum 
slopes ^5 than —2 (as indicated by their J125 — i^ieo colors), with a 
median < —2.5. This suggests that 2; > 8 galaxies are not only dust 
free but also probably have very young ages or low metallicities (see 
also Bouwens et al. 2010). 



binned on a 0.06''-pixel scale. The drizzling was done 
using a modified version of multidrizzle (Koekemoer et 
al. 2002). 4(7 outliers were rejected. 

Given that the WFC3/IR instrument is still relatively 
new, we initially made our own estimates of the photo- 
metric zeropoints by performing PSF-matched photome- 
try on sources present in both the new observations and 
the HUDF NICMOS observations (e.g., Thompson et al. 
2005; Bouwens et al. 2008; Oesch et al. 2009). The ze- 
ropoints derived were consistent (<0.05 mag) with the 
official STScI values, so we elected to use the STScI val- 
ues 26.27, 26.25, and 25.96 mag for the Yios, ^125, and 
i7i6o bands, respectively. The approximate 5cr depths 
of the F, J, and R WFC3/IR images are 28.8, 28.8, 
and 28.8 mag, respectively, in 0.4''-diameter apertures, 
-1.5 mag deeper than the NICMOS data over the HUDF 
(Thompson et al. 2005). These depths were estimated 
by measuring the noise statistics in apertures of various 
size. The FWHM of the PSF in our WFC3/IR near- 
IR images is —0.16''. For reference, the HUDF optical 
^435^606^775^850 data (Beckwith et al. 2006) reached to 
29.4, 29.8, 29.7, and 29.0 AB mag (5cr: 0.35''-diameter 



UDFy-42886345 
I • 




UDFy-43086276 



Fig. 2. — Vboe ^775 2^850 >^105 -^125^^160 image cutouts {2 A" x 2 A" 
on a side) of the five z ~8 Yios-dropout candidates we identified in 
the ultra-deep HUDF WFC3/IR observations. All five candidates 
are detected at 5cr in both the J125 and iifieo bands, but show 
no detection (> 2cr) in the ultra-deep optical HUDF Veoe, '^775, or 
^850 band images. 

apertures) and had PSF FWHMs of -0.10^'. 

3. OBJECT DETECTION AND VERIFICATION 

3.1. Catalog Construction 

Our procedure for doing object detection and pho- 
tometry is identical to that used in previous work (e.g., 
Bouwens et al. 2007; Bouwens et al. 2008) and is per- 
formed using SExtr actor (Bertin & Arnouts 1996: run in 
dual-image mode) on the registered data. Object detec- 
tion is done from the coadded J125 and i^ieo-band image 
(explicitly, using the square root of the image: Sza- 
lay et al. 1999), both of which are redward of the break 
for Yios -dropout galaxies. After smoothing the optical 
data to match the WFC3/IR PSFs, colors are measured 
using Kron (1980)-style photometry in small apertures 
that scale with the size of the source (for Kron factors 
of 1.2). Flux measurements in these small apertures are 
corrected up to total magnitudes using the flux in larger 
scalable apertures (Kron factors of 2.5). This correction 
to total magnitudes (see e.g. Figure 5 from Coe et al. 
2006) is done on a source by source basis, based on the 
square root of image (approximately proportional to 
the coadded flux). Finally, the total magnitudes were 
corrected by 0.1 mag to account for light on the wings of 
the PSF. 

3.2. Dropout Selection 

Sources are selected over 4.7 arcmin^ with the deep- 
est WFC3/IR observations (^10% of this area is not 
at the full depth). The selection criteria we adopt for 
identifying 2; ~ 8 Yios dropout galaxies are simple ana- 
logues of the criteria used to select Lyman Break galax- 
ies at lower redshift (e.g., Steidel et al. 1996; Giavalisco 
et al. 2004; Bouwens et al. 2007). That is, we require 
galaxies to show strong breaks (Y105 — J125 > 0.8) at 
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TABLE 1 

' 8 Yios-DROPOUT CANDIDATES IDENTIFIED OVER OUR ULTRA-DEEP HU DF 
WFC3/IR FIELD. 



Object ID 


R.A. 


Dec 


^160 


^105 - ^125'' 


Jl25 — ^160 


UDFy-42886345 


03:32:42.88 


-27:46:34.5 


28.0±0.1 


1.1±0.4 


0.0±0.2 


UDFy-37796000 


03:32:37.79 


-27:46:00.0 


28.0±0.1 


2.1±0.9 


-0.3±0.2 


UDFy-38135539 


03:32:38.13 


-27:45:53.9 


28.1±0.1 


>2.1 


0.2±0.2 


UDFy-37636015 


03:32:37.63 


-27:46:01.5 


28.4±0.1 


>1.8 


0.1±0.3 


UDFy-43086276 


03:32:43.08 


-27:46:27.6 


29.0±0.2 


1.1±0.5 


-0.5±0.3 



Lower limits on the measured colors are the la limits. Magnitudes are AB 
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Fig. 3. — Surface density of z ~ 8 lios-dropouts found in our 
ultra-deep HUDF WFC3/IR field versus ifieo-band magnitude. 
For comparison, the surface density of Yios-dropouts one would 
expect assuming no evolution from 2; ~ 6 and 2; ~ 7 is also plotted 
(shown with the red and magenta lines, respectively: see §4.1). 
Also plotted (dotted black line) is the expected surface density if 
one extrapolates the z~4-6 LF results of Bouwens et al. (2008) to 
z ~ 8, i.e., with a (j)* of ~0.0011 Mpc~^, a of —1.74, and a fainter 
value of M*. 

the redshifted position of Lya at z ~ 8 and to be blue 
redward of the break, i.e., (J125 — i^ieo < 0-5) and 
(-^125 — ^i6o) < 0.2 + 0.12(yio5 — J125) to exclude intrin- 
sically red galaxies at lower redshift (Figure [1]). We also 
require our Yios-dropout candidates to show no detection 
(< 2a: color-measurement aperture) in all bands blue- 
ward of the dropout band, i.e., ^775^850- Sources 
showing a > 1.5(7 detection in > 1 optical band were 
also eliminated. All of our Yios-dropout candidates were 
required to be 5.5cr detections in the J125 band to en- 
sure they corresponded to real sources. We elected to 
use a 5.5<j criterion to be conservative for these early 
WFC3/IR data. 

The most stringent aspect of the current selection 
is our requirement that sources be undetected in the 
ultra-deep HUDF optical data (which reaches to 31.5- 
32.0 at 1(7 for most sources). From our simulations 
(§3.4), this requirement eliminates almost all contamina- 
tion from z ^ 6 galaxies, and hence we can use a modest 
(^105 — -^125) > 0.8 break to select z ^ 8 galaxies, without 
significant contamination concerns. 

3.3. Yio^-dropout Sample 

We identify five sources that satisfy our F-dropout cri- 
teria (see Table [1]). The position of these sources in the 




-16 



1700.AB 



Fig. 4. — Constraints on the rest-frame UV LF at 2; ~ 8 
from the present Yios -dropout searches over our ultra-deep HUDF 
WFC3/IR field (solid red circles, la error bars and upper limits). 
Also included is our UV LF determination at z ~ 7 from a z- 
dropout search over the same ultra-deep WFC3/IR field (magenta 
lines and points: Oesch et al. 2010a) and also the Bouwens et al. 
(2007) UV LF determinations at z ~4,5 and 6 (shown in blue, 
green, and cyan lines and points). While the other UV LFs shown 
here were derived at slightly different rest-frame wavelengths, their 
Muv,AB at 1700Ais essentially the same (< 0.1 mag: since star- 
forming dropout galaxies are very fiat in fu). The dotted red line 
gives one possible Schechter-like LF that agrees with our search 
results and adopts the Bouwens et al. (2008) scaling of M* with 
redshift. 

^105 - -^125 / -^125 - ^160 two-color plane is given in Fig- 
ure [TJ^ The J125 — HiQo colors of our Yio5-dropout can- 
didates are very blue in general, corresponding to UV- 
continuum slopes of < — 2 (where fx oc A^), with a 
median /3 < —2.5 (though we emphasize the number of 
sources are still small and the inferred /3's may be af- 
fected by uncertainties in the photometry). The observed 
/3's are blue enough that z ^ 7 galaxies must be essen- 
tially dust free, and possibly also have very young ages 
or metallicities (see also Bouwens et al. 2010). 

Cutouts of the sources are provided in Figure [2l The 
candidates have Hiqq^ab band magnitudes of ~28. 0-29.0 
AB mag, within one magnitude of our sensitivity limit. 
Such sources could not have been found to date since 
they are fainter than could be probed with other data 
sets. This illustrates the importance of the very deep 
near-IR data being collected as part of this program. 

All five candidates have apparent half-light radii of 
~0.15'' (~0.7 /i~^kpc) - measured using SExtractor - 
not much larger than the PSF. The 4 candidates - for 
which crowding is not a concern - do not show signifi- 

^ It is encouraging that recently both McLure et al. (2010) and 
Bunker et al. (2010) have independently identified a very similar 
set of z ~ 8 candidates in the first epoch HUDF09 WFC3/IR 
observations over the HUDF. 
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cant (>2cr) detections in the IRAC data over the GOODS 
fields either individuahy or when stacked (Labbe et al. 
2010). The non-detection of the candidates in the IRAC 
data is not particularly surprising given the much shal- 
lower depths of the IRAC data (<27.0 AB mag at 2a) 
relative to the WFC3/IR data. 

3.4. Contamination Corrections 

The only meaningful source of contamination for the 
present sample are sources that enter the selection via 
photometric scatter. A simple estimate of the likely con- 
tamination can be obtained by adding noise to the color 
distribution observed for ~25-26.5 AB mag galaxies. The 
advantage of using this color distribution for the simu- 
lations is that the distribution is realistic (being taken 
from the observations), has a higher S/N than at fainter 
magnitudes, and does not not include any z > 7.5 galax- 
ies. We find ^0.2 contaminants per field, suggesting such 
a source of contamination is minimal (^4%). 

The < 4% contamination estimate implicitly in- 
cludes the contribution of z~ 1.5-2 Balmer-Break Galax- 
ies (BBG) sources scattered into our selection. We ex- 
pect the explicit contribution to be small given the lack 
of faint sources with V — J > 1.5, J — H > 0.15 colors 
expected for BBGs (arbitrary reddenings, metallicities: 
R.J. Bouwens et al. 2010, in prep). 

Other sources of contamination are not important for 
this selection. For example, given that each source is 
detected at >5.5cr in the J125 band and >4cr in the 
HiQo band, no contamination from spurious sources is 
expected. Contamination by SNe is also not important, 
given that the Yios -^125^160 data were taken over the 
same 12-day window. Finally, contamination by stars 
is also unlikely. Not only are T dwarfs rare over the 
CDF-South, with a surface density ^0.04 arcmin"^ (e.g., 
Bouwens et al. 2008) and therefore unlikely to be found 
in our 4.7 arcmin^ field, but also all 5 of our candidates 
appear to be extended (having SExtractor stellarity pa- 
rameters <0.3: see also Oesch et al. 2010b) and hence 
not likely to be stars. 

4. RESULTS AND IMPLICATIONS 
4.1. Expected Numbers 

To interpret the results of the present z ^ 8 F105- 
dropout selection, it is useful to estimate how many Y105- 
dropouts we might have expected if the UV LF showed 
no evolution from z ^ 6 and z ^ 7. We estimate the 
numbers by creating galaxy catalogs according to the 
model LFs, adding artificial galaxies to the data, and 
then processing the images and doing the selection in 
the same way as for the real data. For these simulations, 
we model the pixel-by-pixel profiles of galaxies at z^S 
with similar-luminosity galaxies from the z ^ 4 Bouwens 
et al. (2007) HUDF B dropout sample, but scaled in 
size to match the observed z)~^ size-redshift scaling 
(Oesch et al. 2010b; Bouwens et al. 2004, 2006; Fergu- 
son et al. 2004). Star- forming galaxies at z ~ 7 — 9 are 
assumed to have mean /7V-continuum slopes P of —2.5, 
with a la scatter of 0.4, to match the apparent colors of 
Yios-dropouts in our sample (see Figure [1]). 

Adopting the z ^ 6 z-dropout LF from Bouwens et al. 
(2007: see also McLure et al. 2009) and assuming no evo- 
lution to higher redshift, we predict 20±6 Yios-dropouts 



TABLE 2 

UV Luminosity Densities and Star Formation Rate 
Densities.^ 







\ogiQ SFR density 


Dropout 


(ergs s 


(M odot Mpc-3 yr-1) 


Sample 


< z> Hz-i Mpc-3) 


Uncorrected Corrected^ 


z 


6.8 25.73±0.16^ 


-2.17 ±0.16 -2.17 ±0.16 


Y 


8.2 25.18±0.24 


-2.72 ±0.24 -2.72 ±0.24 



Integrated down to -18.3 AB mag, or 0.08 L^^^. 
^ The dust correction is taken to be 0, given the very blue /3's (see 
also Bouwens et al. 2009; Bouwens et al. 2010) 

^ Using the z 7 UV LF derived from the same ultra-deep 
WFC3/IR field as this study (Oesch et al. 2010a) 
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Fig. 5. — Determinations of the UV luminosity density and SFR 
density, integrated to 0.08 L*^^ (—18.3 AB mag) as appropriate 
for the z 8 sample, as a function of redshift. The large red 
circle shows the constraints we can set on this density at z ~ 8 
from the current lios -dropout search (see also §4.3). The lower 
set of points (and blue region) show the SFR density determintion 
inferred directly from the UV light, and the upper set of points 
(and orange region) show what one would infer using dust correc- 
tions inferred from the f/y-continuum slope measurements (e.g., 
Bouwens et al. 2009; Bouwens et al. 2010). Note that the dust 
correction is essentially zero at z > 6. Also included on this fig- 
ure are the determinations at 2; ~ 7 from the HUDF WFC3/IR 
z-dropout search (Oesch el. 2010a: solid red square), the Bouwens 
et al. (2007) determination at z ~ 4 — 6 (open red circles), the 
Reddy &; Steidel (2009) determinations at z ~2-3 (green crosses), 
and the Schiminovich et al. (2005) determinations at z < 2 (black 
hexagons). A systematic increase in the SFR density from z ~ 8 
to z ~ 2 is clear. 

in our WFC3/IR field. The z ^ 6 predictions are shown 
as a function of i^ieo-band magnitude in Figure [3] and 
compared with the observations. We expect ~70% un- 
certainties in these numbers due to small number statis- 
tics and large-scale structure variations (assuming a 4.7 
arcmin^ survey field, Az ^ 0.8 redshift selection window, 
and ~ 10~^ Mpc~^ source density: e.g., Trenti & Sti- 
avelli 2008). The observed numbers are > 2a lower than 
expected from the well-determined z ^ 6 LF. Predicted 
surface densities can also be made using the Bouwens et 
al. (2008) z ^ 7 LF (also shown), but the uncertainties 
are large given the small sample in that paper (and the 
new WFC3/IR observations permit much better z ^ 7 
LF estimates). 

4.2. Implications for the z ~ 8 UV LF 

The fact that the observed numbers are lower than pre- 
dicted assuming no-evolution from z ^ Q suggests that 
the UV LF continues to evolve from z ~ 8. So, it is in- 
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teresting to extrapolate the LF results from Bouwens et 
al. (2008) to z - 8 - giving M^y = -19.45, ^* - 0.0011 
Mpc~^, and a ~ —1.1 A - and see what we find. Per- 
forming this exercise, we predict that 5 Yios-dropouts 
would be found in the present search {shown with the 
dotted black line in the lower panel of Figure |4]) . This is 
in good agreement with the observed results. 

We can further quantify the overall magnitude of this 
evolution. Using the Bouwens et al. (2008) LF param- 
eterization as a guide, we fix a = —1.74, (/)* = 0.0011 
Mpc~^, and then derive confidence intervals on M^y. 
For our Yios-dropout search, we estimate that M^y = 
— 19.5 ± 0.3 AB mag. This is significantly fainter than 
the Mljy = -20.2 ± 0.2 AB mag estimated at z 6 or 
the Mljy -21 AB mag at z 4. 

Of course, given the size of the sample and lack of com- 
plementary wide- area searches for bright z ^ S sources, 
it is difficult to constrain the shape of the 2; ~ 8 UV LF. 
Therefore, we simply consider the stepwise LF at z ~ 8. 
We divide our dropout sample into 0.5 mag bins, com- 
pute the equivalent absolute magnitudes in each of these 
bins, and then divide the observed number of sources 
in each bin by the effective selection volume, which are 
estimated using the same simulations described in §4.1. 
These stepwise LF determinations are presented in Fig- 
ureSl with the LFs at (Bouwens et al. 2007; Oesch 

et al. 2010a) shown for context. The la upper limits on 
the volume density of luminous z ^ S sources are also 
shown. It would appear that the UV LF only shows very 
weak evolution at low luminosities (~ —18.3 AB mag). 
This is in contrast to the dramatic evolution observed at 
the bright end from z ~ 7 to z ~ 4 (see e.g. discussion 
in Shimasaku et al. 2005; Bouwens et al. 2008). 

4.3. Constraints on the UV Luminosity Density/SFR 
density at z ^ 8 



Finally, we calculate the luminosity densities (and un- 
obscured SFR densities) at z ~ 8 implied by these con- 
straints on the rest-frame UV LF. For the luminosity 
density at z ~ 8, we simply integrate the stepwise z ^ S 
LF shown in Figure [H We convert these UV luminosity 
densities into the equivalent unobscured SFR densities 
using the Madau et al. (1998) prescription. The results 
are presented in Figure [5] and Table [2l Also presented 
are the star formation rate densities inferred. The dust 
correction is taken to be 0, given the very blue /3's (see 
also Bouwens et al. 2009; Bouwens et al. 2010). 

As these results demonstrate, the remarkable improve- 
ment in the sensitivity and the "discovery efficiency" 
(area gain x sensitivity gain) of WFC3/IR has enabled 
HST to cross a threshhold. HST can now find star- 
forming galaxies at z~8-8.5. The existence of such galax- 
ies and the active star formation implied at even higher 
redshifts z > 10 provides a striking framework for future 
detailed JWST observations. 
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